Collective Behavior of the Closed-Shell Fermi Gas 
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We propose an unconventional description for the ground state and collective oscillations of the 
two-component normal Fermi gas with two-body zero-range interactions. The many-body problem 
can be accurately reduced to a linear, one-dimensional Schrodinger equation in a single collective 
coordinate, the hyperradius R of the iV-atom system which is the root mean square radius. The 
calculated properties of the Fermi gas ground state are shown to agree accurately with results from 
the Hartree-Fock (HF) approximation over a wide range of interspecies scattering lengths. The 
breathing mode excitation frequency deviates qualitatively from HF predictions, but we show that 
this reflects a failure of the HF approach for this observable. 



Tremendous excitement has ensued in the ultracold 
physics community since dilute 2-component Fermi gases 
have recently been shown to exhibit BCS-like pairing 
physics. |l| This pairing occurs at sufficiently low tem- 
perature T << Tf, sufficiently attractive interatomic 
scattering length a, and sufficiently high density p = 
k F /(3ir 2 ), where kp is the Fermi wavenumber. This in- 
triguing physical regime of the BEC-BCS crossover arises 
in the range kpa < — 1. The DFG currently studied 
in this region experimentally are in a metastable state 
whose decay to a more compact geometry in a different 
region of configuration space is inhibited by a multiparti- 
cle tunnelling barrier. Comparatively little attention has 
been given to the behavior of the low-temperature normal 
Fermi fluid with attractive interactions. In this Letter 
we discuss the behavior of this normal Fermi fluid from 
a viewpoint that departs from conventional mean-field 
theory. This theoretical framework and subsequent im- 
provements may prove useful for analyzing other many- 
body systems, because in contrast to the familiar mean- 
field theories that require the solution of inherently non- 
linear differential or integro-differential equations, the 
present treatment utilizes only linear Schrodinger equa- 
tions throughout. 

Our treatment is related to the hyperspherical coordi- 
nate treatments that have been applied to finite nuclei, 
as in the "K-harmonic" method reviewed by Smirnov 
and Shitikova. [2j The formulation is a rigorous varia- 
tional many-body calculation, aside from the usual lim- 
itations of the assumed pairwise Hamiltonian with zero- 
range interactions. As was found in a similar study 
for bosons, the resulting hyperspherical coordinate de- 
scription again provides a linear, one-dimensional, effec- 
tive Schrodinger equation in the hyperradius. In this 
initial study of the usefulness of this methodology, we 
restrict ourselves to the study of fermionic systems with 
completely filled shells. This is done for analytical and 
computational simplification; the same approach should 
be equally applicable to open-shell fermi systems, with 
modest extensions. 

Consider N identical fermions confined either magnet- 
ically or optically in a spherically symmetric harmonic 



oscillator trap with angular frequency to. Half of the 
fermions are in one spin substate, while the other half 
are in a second spin substate. The full TV-body time- 
independent Schrodinger equation for this system reads 
Hip = Eip, with 
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Here m is the mass of the fermions under consideration 
and Ui n t is the pairwise atomic interaction potential. 
We adopt the Fermi pseudopotential Ui n t ~ ?j) = 
4:Trh 2 a , 
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where a is the two body s wave scatter- 



ing length. Next we transform to a set of hyperspherical 
coordinates. The hyperradius for this system of equal 
mass particles is the root mean square distance of the 

1 /2 

atoms from the trap center, R = (y2r 2 /N) . The re- 
maining 3N — 1 dimensions are described by a set of an- 
gular coordinates collectively denoted by fi. In addition 
to the ordinary set of 2N spherical polar coordinates for 
each particle {{0i,4>i)}^ =1 , we follow Ref. [3] in defining 
N — 1 hyperangles a.;, £ = 1,2, 3, ...,N — 1. as: 
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tancv, ee (^Vf) 1/2 /r i+ i. 
In this coordinate system the kinetic energy becomes 
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where M — Nm is the total mass of the system and 
the squared grand angular momentum operator A 2 is 
defined in the standard fashion 0, • The trap potential 

is purely hyperradial, i.e. J^^i -^m^r 2 — -Muj 2 R 2 . 

Multiplying ip by i?( 3Ar_1 )/ 2 allows us to eliminate first 
derivative terms. One can view the rest of the calculation 
as the simplest approximation to Macek's adiabatic hy- 
perspherical description [5]or as a one-term truncation of 
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an expansion of the exact ip (R, O) into the complete, or- 
thonormal set of eigenfunctions of the operator A 2 , mul- 
tiplied by unknown radial functions that are optimized 
variationally. These eigenfunctions, or "hyperspherical 
harmonics" (HHs, see 00) can be expressed analyti- 
cally in terms of Jacobi polynomials, for any number of 
particles. Their eigenvalue equation is 

A 2 Y Xu (ft) = A (A + 3N - 2) Y Xv (O) 

where A = 0,1,2,.. is the order of the harmonic. The 
second index, v, stands for the further set of N — 2 quan- 
tum numbers required to distinguish the (usually quite 
large) set of degenerate states having any chosen value 
of A. This allows for the choice of a trial wave func- 
tion if) (R, fi) = G (R) $a (to) where $ A (to) is a linear 
combination of HHs (multiplied by appropriate spinors) 
with grand angular momentum A, which satisfies the an- 
tisymmetry requirements of this fermionic system. Gen- 
erally the coefficients of this expansion are very diffi- 
cult to obtain; in previous studies they have usually 
been calculated recursively 0, IE using a coefficient of 
fractional parentage (cfp) expansion. To eliminate the 
need for a demanding cfp calculation with hundreds or 
thousands of particles, we restrict our treatment here 
to only magic numbers of particles, i.e., to filled energy 
shells of the noninteracting spherical harmonic oscilla- 
tor. Because the TV-particle oscillator has a separable 
Schrodinger equation in either independent particle co- 
ordinates or hyperspherical coordinates, any nondegen- 
erate, normalized, antisymmetric iV-body state must be 
identical in the two coordinate systems. The nondegener- 
ate wavefunction for n completely filled shells (n = 1,2...) 
can thus be written as a Slater determinant of the in- 
dependent particle states, which allows us to construct 
$a (0) as: 
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Here R ni ii (>t) is the radial solution to the independent 
particle harmonic oscillator for the ith particle given 



(r/iy 



by rR nl (r) = N nl exp (-r 2 /2Z 2 ) (r/l) t+1 L l ,+ 1/2 

where (r) is a modified Laguerre polynomial with 
/ = \Jhjmuj. yi imi (w) is an ordinary 3D spherical 
harmonic with LUi as the spatial solid angle for the 
ith particle, \m Si ) is a spin ket that allows for two 
spin species of atoms, |f)and ||). R ( - 3N ~ 1 ^ 2 F X (R) = 
A A exp (-i? 2 /2L 2 ) (i?/£) A+3(JV - 1)/2+1 is the nodeless 
hyperradial solution for N noninteracting particles in an 
oscillator with hyperangular momentum A where A x a 
normalization factor and L = l/\/N. With this trial 
wavefunction, the expectation value of the Fermi pseu- 
dopotential yields an i?-dependent effective potential for 
the function G(R). From known properties of the S func- 
tion and the standard properties of matrix elements of 



two-body operators between two Slater determinantal 
wave functions, the desired diagonal potential matrix el- 
ements can be shown to have the form (in harmonic os- 
cillator units for both energy and radius) 
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The constant C\ only depends on the number of atoms 
in the system, and is found by examining the diagonal 
matrix element of this hyperradial potential in the state 
F\ (R). Integration over the hyperradius allows us to 
solve for C x : 



_ N (N — l)l 5 T(K + 3/2) 
~ 2T (K) 

x (F A (R) $ a (0)| 5 (n - r 2 ) \F\ (R) $ A (0)) 



(3) 



Here K = X + 3 (N - 1) /2 and F x (R)<P x (n) 
is the fully anti-symmetrized, non-interacting 
atom wave function. Thus the matrix element 
(F x (R)^ x (Q,)\S(f 1 -f 2 )\F x (R)^ x (Q)) is a sum of 
integrals over products of Laguerre polynomials in 
independent particle coordinates. Every diagonal ma- 
trix element coefficient C x for the n-th complete filled 
shell can be written as a rational number multiplied 
by \/2/7r 2 . For instance, the n — 2 filled shell state 
has N — 8 particles, and in its ground state, grand 
angular momentum A = 6, and C x has the value 
(134217728/570285)v^/tt 2 . For the n = 3 filled shell 
with 20 particles, C x involves the ratio of two integers, 
the larger of which is in the numerator with 34 digits. 
In view of the rapid increase as the number of particles 
is increased, we focus on the asymptotics in the large 
TV limit. The following functional form is found to be 
reasonably accurate for all N beyond about N = 20 
particles, namely C X (N) -► 0.16802iV 7 /2y2/7r 2 . The 
fractional error in this approximation is smaller than 
1% for N = 20 particles, and smaller than 0.1% for 
the 30th filled shell with 9920 particles, and beyond. 
The computational strategy has enabled us to carry 
out explicit calculations for as high as the 50th magic 
number, which is a degenerate fermi gas with 44,200 
atoms. This is high enough to extract the leading 
behavior in the high- TV limit. 

This effective potential now gives a linear, one- 
dimensional Schrodinger equation in the hyperradius 
h [R^ N - 2 )/ 2 G (R)] = E [R^ N ' 2 y 2 G (R)\ where h (in 
oscillator units) is given by 
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It is informative to recast this radial Schrodinger equa- 
tion into dimensionless form with E = EnjE' and R = 
^/{R 2 ) NI R', where E NI = (A + 3JV/2) and {R 2 ) NI = 
(3/2 + X/N) are the energy and expectation values of 
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R 2 of the non-interacting ground state (in oscillator 
units). The resulting linear effective one-dimensional 
Schrodinger equation is given by 



(5) 

where m* = (A + 3A/2) Z , and, in the large A limit, 

4ttciC\ 



V e f f (R) is given by 

1 I 



Veff (#) 



2R' 2 2 
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I (A + 3A/2) 5/2 i?' 3 



The value of A for each filled shell ground state can 
be ascertained by comparing the energy of A nonin- 
teracting particles in a harmonic oscillator potential, 
written in each coordinate system. In the indepen- 
dent atom representation the energy is given by E = 

huj Y,f=i {h + 2n,) + 3A/2 . The energy of A nonin- 
teracting particles, described using hyperspherical quan- 
tum numbers, with grand angular momentum A and 
X = 0,1,2,... hyperradial nodes in F\(R), is given by 

E = hu ^A + 2\ + ~jT^J ■ The ground state has \ = 0, 

while the first breathing mode excited state has x = 1. 
After we equate the two expressions for the ground state 
energy, we obtain A = Yli=i (h + 2ri»). In terms of 
the number of complete filled shells, n, namely N(n) = 
n(ra + l)(n + 2)/3, and A(n) = (n - l)n(n + l)(n + 2)/4. 
To a good approximation at all N, A(A) ~ (3A r ) 4/3 /4 - 
3A/2 + A^ 2 / 3 /192 1 / 3 , which becomes exact in the limit 
N — > oo. Similarly, the Fermi wave number, kf, may be 
shown to be k f = ^2{n+l/2)/l -> y/2 (3 1 / 6 ) N^ e /l. 
This all shows that in the large- N limit 
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where C ^ 0.16802 (64Vl2/8l7r) . Note that the only pa- 
rameter in this asymptotic large-TV potential curve is 
in fact fc/a. At N — > oo we see that m* — > A 8 / 3 , 
whereby the second derivative term in equation 0] is neg- 
ligible. This means the squared zero-order ground state 
wave function in the large A limit is approximated by 
[R>m-i)/2G (R')] 2 w SiR'-R'^) where R' min is the 
location of the minimum of V e ff (R 1 ), if it exists. This 
recognition of the strong hyperradial localization for in- 
finite dimensional systems is also used extensively in di- 
mensional perturbation theory. Ground state expecta- 
tion values of quantities involving the hyperradius can 
be extracted by simply evaluating them at R' min (e.g. 
(i?' 2 > = R'X and E = V e f f (R'^)). 

The effective potential in equation^ is similar in form 
to the analogous hyperradial potential obtained for the 
Bose-Einstein condensate (BEC) in Ref. 0. The main 
differences between these two potentials stem from the 
antisymmetry requirements; whereas the constant in the 



interaction potential acted as N 3 / 2 for a BEC, we 
find for the DFG an A 7 / 2 behavior for A >> 1. On the 
other hand, the coefficient of the R~ 2 effective centrifu- 
gal repulsion increases as A 5 / 3 , as opposed to N for the 
BEC case, which is a manifestation of the Fermi pressure 
in the DFG. The hyperspherical results for the ground 
state energy and expectation value of R' 2 are very close 
to the HF predictions, for both repulsive and attractive 
interactions, as is evident from Fig. ^ 
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FIG. 1: The ground state energy (solid line) and (-R 2 ) (dashed 
line) calculated using the hyperspherical method are com- 
pared with the Hartree-Fock results (circles and squares, re- 
spectively) , as functions of fe/a. All calculations on this figure 
were carried out for 240 atoms. 

The effective potential energy V e ff for values of kta > 
0, i.e. for repulsive interactions a > 0, shows only one 
global minimum in V e ff {R'), & n d the condensate energy 
increases as does the overall size of the DFG. For val- 
ues of kf a < 0, where the interspecies interaction is at- 
tractive, the energy of the ground state decreases below 
the non-interacting energy as the atoms are pulled in 
together by the attractive interaction and by the trap. 
Provided kfa is not too negative, the repulsive 1/i? 2 
kinetic energy term is sufficiently large and positive to 
stabilize the gas in a local hyperradial potential mini- 
mum. But when kf a decreases to sufficiently negative 
values, the local minimum in V e ff {R) disappears alto- 
gether at a critical value we denote a c , given in the large 
iV limit by fc/a c ~ —1.21. This prediction is in good 
agreement with the value of the critical scattering length 
found by the authors of Ref. [10( beyond which there 
was no self-consistent Hartree-Fock solution. This value 
for the critical scattering length differs by less than 8% 
from the value of aN 1 / 6 jl = —.66 found by the authors 
of 11]. Thus we see that there are two regions of interest 
kfd c < kfa < and kfa < kja c . 

Figure shows Veff (R r ) for five equally spaced values 
of kfa about fc/a c ranging over fc/a c — 0.5 to fc/a c + 0.5, 
in steps of 0.25. The behavior of V e ff here resem- 
bles the effective hyperradial potential for N identical 



bosons. BECs formed with a < a c . were predicted to 
collapse by a number of studies 0, 0, , and the cor- 
responding phenomenon of the bosenova was observed 
experimentally, [lil Il5| 



V ef f 




FIG. 2: The effective potential for five values of fc/a spaced 
evenly from fc/a c — 0.5 to fc/a c + 0.5. The bold curve is at 
kfa = fc/a c , and the dotted curve is the non-interacting po- 
tential. 

Moreover, as fc/a is lowered, the local minimum of 
V e ff (R 1 ) becomes "softer", as the curvature at the local 
minimum decreases. This suggests that the lowest exci- 
tation frequency (the frequency of the breathing mode) 
should decrease as a is made increasingly negative. Fig- 
urc[3]compares the excitation frequency calculated in our 
hyperspherical method with the lowest one predicted us- 
ing the sum rule method. |l6| in which the desired mo- 
ments were calculated using the Hartree-Fock ground 
state. The calculated HF frequencies for the lowest 8 ex- 
citations caused by the same R 2 breathing mode operator 
are also shown. The hyperspherical frequencies are in ex- 
cellent agreement with the sum rule predictions, but the 
Hartree-Fock frequencies are qualitatively different and 
apparently erroneous. Our hyperspherical description 
allows for the calculation of collective monopole excita- 
tions in which every particle in the cloud oscillates in and 
out at the characteristic frequency, while Hartree-Fock 
predicts only single-particle excitations in which just one 
orbital of the DFG gets excited. 

In summary, we have developed a fully antisym- 
metrized, variational description of N fermionic atoms, 
with zero-range interactions in a spherical harmonic 
trap. Truncation to the lowest antisymmetric hy- 
perspherical harmonic reduces the problem to a linear, 
one-dimensional Schrodinger equation in the hyperradius 
with an effective potential. In the large atom num- 
ber limit, it only depends on the quantity fc/a. The 
ground state energy and radius predicted from this lin- 
ear quantum mechanical treatment map accurately onto 
those predicted with the Hartree-Fock method, while the 
breathing mode frequencies agree with the sum rule and 
are a qualitative improvement over the calculated HF 
frequencies. This picture predicts a ground state energy 
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FIG. 3: The lowest excitation energy is shown as a function 
of fc/a. The negatively sloping dashed lines are the lowest 8 
frequencies predicted by Hartree-Fock. The positively sloping 
solid line is the hyperspherical prediction and the circles are 
the predictions from the sum rule calculation. All of these 
calculations were carried out for the closed shell system having 
a total of 240 atoms. 



that collapses to — oo at fc/a < —1.21, in a manner iden- 
tical to the physics of the "bosenova". 0But the full in- 
terrelationship between this phenomenon and the physics 
of pairing in the BEC-BCS crossover region goes beyond 
the scope of the present study and will be explored in a 
subsequent publication. 
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